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We present a symmetry-based analysis of competition between different gapped states that have 
been proposed in bilayer graphene (BLG), which are all degenerate on a mean field level. We classify 
the states in terms of a hidden SU(4) symmetry, and distinguish symmetry protected degeneracies 
from accidental degeneracies. One of the states, which spontaneously breaks discrete time reversal 
symmetry but no continuous symmetry, is identified as a Quantum Anomalous Hall (QAH) state, 
which exhibits quantum Hall effect at zero magnetic field. We investigate the lifting of the accidental 
degeneracies by thermal and zero point fluctuations, taking account of the modes softened under 
a renormalisation group procedure (RG). Working in a 'saddle point plus quadratic fluctuations' 
approximation, we identify two types of RG-soft modes which have competing effects. Zero point 
fluctuations, dominated by 'transverse' modes which are unique to BLG, favor the QAH state. Ther- 
mal fluctuations, dominated by 'longitudinal' modes, favor a SU(4) symmetry breaking multiplet of 
states. We discuss the phenomenology and experimental signatures of the QAH state in BLG, and 
also propose a way to induce the QAH state using weak external magnetic fields. 

PACS numbers: 



I. INTRODUCTION 

The Quantum Anomalous Hall (QAH) insulator is a 
state of matter where spontaneous breaking of time re- 
versal symmetry produces (integer) quantum Hall effect 
in the absence of any external magnetic field. First pre- 
dicted in 1988 1 , the QAH state has never yet been ob- 
served. In the recent literature on interaction driven 
topological insulators^, the elusiveness of the QAH state 
has been ascribed to fluctuations, which typically disfa- 
vor the QAH state with respect to a Quantum Spin Hall 
(QSH) state, which is degenerate with the QAH state on 
a mean field level. Here we point out that the fluctua- 
tions which govern the competition of different gapped 
phases proposed in bilayer graphene (BLG)^~— are domi- 
nated by the modes not present in the models^ 3 -, leaving 
open the door to formation of a QAH state at zero field 
in BLG. Also, we will propose a mechanism for inducing 
the QAH state using external fields. 

The theoretical literature on BLG predicts instabili- 
ties to numerous strongly correlated states, which are 
gapped 4 ^— or gaplessi depending on the way the electron- 
electron interaction is modeled. The numerous gapped 
states predicted in the literature are all degenerate at the 
level of mean field theory 4 ^, and have the same instabil- 
ity threshold under one loop RG-. The relation between 
these different states, and their experimental signatures, 
have not yet been understood. 

Meanwhile, recent experiments indicate that the 
gapped state observed in charge neutral BLG in quan- 
tizing magnetic fields - persists down to low fields, cross- 
ing over to another gapped state at zero field^. However, 
the nature of the gapped state at zero field is unknown. 
Hence, clarifying the relation between different gapped 
states and understanding their physical properties is an 
interesting and timely task. 

Here we present a unifying symmetry based analysis of 
strongly correlated states in BLG. The states predicted 



in Refs.[3-[ll are classified according to a hidden SU(4) 
flavor symmetry into symmetry breaking multiplets and 
an SU(4) invariant singlet. The SU(4) singlet is a QAH 
state. The degeneracy of the multiplets and the singlet 
is an artefact of the approximations made in the analy- 
sis, and will be lifted upon taking fluctuation effects into 
account. 

Our analysis of fluctuations in BLG focuses on the ef- 
fect of the modes softened under RG. Those include the 
'longitudinal' fluctuation modes (L-modes) analogous to 
those discussed in Ref-HQ, and also 'transverse' fluctu- 
ation modes (T-modes) which are unique to BLG. We 
find that these two types of modes have competing ef- 
fects: while the L-modes favor the symmetry breaking 
multiplets, the T-modes favor the SU(4) invariant QAH 
state. The zero-point fluctuations are dominated by the 
T-modes, and hence appear to favor a QAH state at zero 
temperature. Meanwhile, thermal fluctuations are domi- 
nated by the L-modes, and favor the symmetry breaking 
multiplets. We speculate that thermal fluctuations may 
drive a phase transition from the QAH state at low tem- 
peratures to a SU(4) symmetry breaking state at higher 
temperatures, and estimate the transition temperature. 
We also discuss the phenomenology of the QAH state, its 
possible experimental signatures, and propose a way to 
further stabilize it using external magnetic fields. 



II. SU(4) SYMMETRY 

In this section we show that within the often-used ap- 
proximation where the difference between interlayer and 
intra-layer interactions is neglect ed 4 -"— i 1 3 ' 1 4 , the interact- 
ing Hamiltonian is invariant under rotations in a suitably 
defined four-dimensional flavor subspace. Specifically, we 
perform a unitary transformation by exchanging the sub- 
lattices A and B in one of the valleys, upon which the 
single particle Hamiltonian becomes identical for all spin 
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and valley species, while the layer and sublattice blind 
interactions are left unchanged. 

Before entering the discussion of the SU(4) invariance 
in BLG, we recall that electronic states in BLG at low en- 
ergy are described by wave-functions on the A and B sub- 
lattices of the upper and lower layers^, and are fourfold 
degenerate in spin and valley. To analyze the structure 
of the Hamiltonian, it will be convenient to combine the 
spin and valley components in a single eight-component 
wavefunction ipa,s,v ( x )j where a is the sublattice (layer) 
index. We shall use the Pauli matrices in sublattice, spin 
and valley space, denoted below by Tj, <Ji and r]i, respec- 
tively. The low energy non-interacting Hamiltonian may 
then be written as 
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where t± = t\ ± iT2- Here m = 0.05m e is the effective 
mass. Because of the presence of rj 3 in Eq.([T]), the sin- 
gle particle Hamiltonian is not invariant under rotations 
of valley components. To bring it to an SU(4) invariant 
form, we perform a unitary transformation on all opera- 
tors 
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This transformation does not act on the spin space, how- 
ever it mixes the layer and valley indices of the wavefunc- 
tion ^a,s,i)( x ) by interchanging the r pseudospin compo- 
nent (layer) in one of the valleys. As a result, r+ and 
r_ are interchanged and T3 changes sign in the 773 = —1 
valley, after which the free-particle Hamiltonian, Eq.([T]), 
becomes identical in both valleys. 

Defining p± = p x ± ip y , the transformed non- 
interacting Hamiltonian takes the compact form 
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where f + and f_ are obtained by transforming r + and t_ 
according to Eq.Q. This single particle Hamiltonian is 
manifestly invariant under SU (4) rotations in the trans- 
formed spin/valley flavor space. 

Meanwhile, electron interactions can be described by 
a many-body Hamiltonian written in terms of p q = 
S P ^p^p+q (the density summed over layers) and A q = 
S p '0p'?3'73'0p+q (the density difference between layers). 
The interacting Hamiltonian, which incorporates a differ- 
ence between interlayer and intralayer interaction^, can 
be written as 

H = J] ^ p H ^ P + ^Y, V + W-q + V- \ X -^ ( 4 ) 



where V + {q) = 2ite 2 / nq is the Coulomb interaction, and 
V- = ire 2 d/ k accounts for the layer polarization energy 
(here d = 3. 5 A is the BLG layer separation). The pp 



term, which is isotropic in flavor space and thus is SU(4) 
invariant, dominates because d is small compared to 



a = H 2 K/me 2 — 10k A, 



(5) 



the characteristic lengthscale set by interactions 5 . Wc 
therefore approximate by neglecting V_, an approxima- 
tion that becomes exact in the weak coupling limit, where 
d/aa — > 0. Under this approximation, the Hamiltonian 
is invariant under SU(4) flavor rotations, generated by 
the operators fji and di. We will henceforth drop the ~ 
symbols for notational convenience, and will refer to the 
operators f , fj and a as r, 77 and a respectively. All op- 
erators are assumed to be transformed operators unless 
specified otherwise. 



III. CLASSIFICATION OF STATES AND 
TOPOLOGICAL PROPERTIES 

In the transformed basis, the mean field Hamiltonian 
for the gapped states described irt^~— may be written as 
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where m = 0.05m e is the effective mass. Here the Pauli 
matrices Tj act on the transformed sublattice space, and 
Q is a 4 x 4 hermitian matrix in the transformed spin- 
valley space (flavor space), satisfying Q 2 = 1. 

Since unitary hermitian matrices have eigenvalues ±1, 
all gapped states can be classified as (M>, M<), where 
M > and M < are the numbers of +1 and —1 eigenval- 
ues of Q respectively. There are three general types of 
states: (2, 2), (3, 1), and (4, 0). There is an additional Z 2 
symmetry associated with the overall sign of Q which is 
absorbed into the sign of A. Following Refs.fHfl5l|. the 
Hall conductance of a state (M>, M<) can be written as 



= (M>-M<)- 



(7) 



where we took into account an additional factor of 2 due 
to the Berry phase 2tt in BLG—. The (4,0) and (3,1) 
states, which have M> =^ M<, thus exhibit a quantized 
Hall conductance at zero magnetic field-the hallmark of 
a QAH state. Because these states have a xy ^ 0, they 
must spontaneously break time reversal symmetry. We 
will henceforth focus on comparing the (4,0) and (2,2) 
states, since the (3,1) states are intermediate between the 
two. We will refer to the (4,0) state as the QAH state, 
but it should be remembered that the (3,1) states are also 
QAH states. In contrast, the (2,2) states have a xy = 0, 
and preserve time reversal symmetry, but instead exhibit 
quantum flavor Hall effect. If we parameterize the flavor 
space by Pauli matrices 77.; and di in transformed valley 
and spin space respectively, then the Q = 03 state is a 
QSH state (Ref.[16j), while the Q =J73 state is a Quan- 
tum Valley Hall (QVH) state (Ref.[l7|). 



3 



(a) QAH state (b) QAH state (c) QSH state 



BA<0 v=4 BA>0 



■ 

c 
LU 



+|z\j 



-4 



4-+ 



« = 3 
w = 2 

M=l 

n = 0' 

H = 

n = -l 

w = -2 
w = ~3 



FIG. 1: Landau level spectrum of the QAH and QSH states. 
Note an anomalous Landau level in the QAH state that has 
no particle- hole-symmetric counterpart. Occupation of this 
anomalous Landau level allows the QAH state (a) to lower its 
energy relative to the states (b,c) at filling factor v = A. 

These states are analogs of the 'topological Mott insu- 
lators' discussed in Refs.[lQ, and as such host topolog- 
ically protected edge states. The counter-propagating 
valley modes for the QVH state were worked out 
in Ref.[17], the co-propagating charge modes and the 
counter-propagating spin modes for the QAH and QSH 
states follow similarly. The protection of edge modes 
is strongest for the (4,0) state due to the unidirec- 
tional, chiral character of these modes. The counter- 
propagating spin currents in the QSH state are protected 
in the absence of spin- flip scattering, while the counter- 
propagating valley currents in the QVH state are pro- 
tected in the absence of intervalley scattering (e.g. by 
short range disorder). 

We note that the above classification of states superfi- 
cially resembles that arising in an entirely different prob- 
lem, namely the Quantum Hall Ferromagnct (QHF) in 
graphene in quantizing magnetic fielcU^. In the latter 
case, however, the integers M > and M< are fixed by the 
electron density, i.e. by filling of the four-fold degenerate 
zeroth Landau level. In the QHF problem, spontaneous 
time-reversal symmetry breaking cannot occur: the ana- 
log of the (4,0) QAH state is a fully filled zeroth Landau 
level exhibiting quantized Hall conductance 2e 2 /h. Fur- 
thermore, in the QHF problem there is no competition 
between states with different M > and M < values, which 
is the main question of interest for us here. 



IV. LIFTING ACCIDENTAL DEGENERACIES 
USING EXTERNAL FIELDS 



transformations, and are hence degenerate within the ap- 
proximations leading to SU(4) symmetry. 

In contrast, the degeneracy of the different manifolds 
(M>,M<) is purely accidental, and may be lifted in the 
presence of a weak SU(4) invariant perturbation. As 
an example, we consider application of a weak trans- 
verse magnetic field B. Incorporated in the Hamilto- 
nian ([5]) through the replacement p — >• p — eA, it pre- 
serves the SU(4) symmetry, and causes the spectrum 
to split into Landau levels^ with an energy spacing of 
order Huj c , where uj c = eB/mc. The Zeeman energy 
2[i,bBgz is not 5/7(4) invariant, but may be neglected 
since tiui c 2/jb_B. When A = and Zeeman terms 
are neglected, the Landau level spectrum is particle-hole 
symmetric and is fourfold degenerate in flavors. 

Crucially, the T non-invariance of the mass term AT3 
means that the Landau level spectrum for the (4,0) state 
is not invariant under B — s> —B and is not particle-hole 
symmetric^. In particular, the zeroth Landau level, 
which has an additional two-fold orbital degeneracy^, 
forms at energy AsignB only, and has no counterpart at 
— Asignf? (see FigOJi,b). This breaking of particle hole 
symmetry can be exploited to induce the (4,0) state using 
magnetic fields. 

We illsutrate this by comparing the energies of the 
(4,0) QAH state and the (2,2) QSH state in external 
magnetic field at filling factor v — 4. In a magnetic 
field, these states are no longer degenerate, because of 
the anomalous Landau level. It is clear from FigQ]that 
the QAH state with the appropriate sign of A (such that 
BA < 0) is favored over the QSH state: 
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where $o is the flux quantum and Fm > .m < is the free en- 
ergy per unit area for a state (M>, M<). This mechanism 
for lifting the degeneracy between QAH and QSH states 
in favor of the QAH state applies to all systems where 
there is such a degeneracy, including the models studied 
in Refs.011]- Of course, at finite B, there is no time re- 
versal symmetry, so the state realized is not a true QAH 
state, but rather is a state showing quantum Hall effect 
at anomalously low magnetic fields, which is smoothly 
connected to a QAH state at B = 0. 

The analysis above is valid only for sufficiently small 
B, when BLG at v = 4 is not far from charge neutrality. 
This is because the excitonic instability that generates 
the gap A (Refs.0-(1]) is suppressed by detuning away 
from charge neutrality. 



The SU(4) symmetry of the Hamiltonian guarantees 
the degeneracy of all states within a given manifold (M> , 
M<), even when the states involved have very different 
physical properties. For example, the QVH state is a 
ferroelectric state which polarizes the layers by charge, 
while the QSH state polarizes the layers by spin and val- 
ley. Nonetheless, the two states are related by SU(4) 



V. SADDLE POINT ANALYSIS 

We now investigate the energy splitting between the 
different manifolds at B — by going beyond a mean field 
approximation, and including the effect of fluctuations. 
We consider BLG in the presence of screened Coulomb in- 
teractions between electrons. A static screening approx- 



4 



imation, ignoring the effects of dynamical screening^, is 
sufficient to understand the main features. In this ap- 
proximation, the interaction is short range, and we can 
write the partition function as a functional field integral 
in Euclidean time, 



Z = J D^Dipexpf- J i 



dxC[ip* (x), ip{x) 



(9) 



where x — (t,r), dx — dtd 2 r, the ip fields are fermionic 
fields, with the Lagrangian 

C = ^(d t + H )rP+ Yl ^HM- (10) 

j,fc=1...8 

Here j, k are combined sublattice and flavor indices, and 
H is the non- interacting Hamiltonian (given by Eq. © at 
A = 0). The coupling constant A represents the statically 
screened Coulomb interaction, which in the RPA model 
takes valued A = l/(4^ohi4), where vq = m/2ir is the 
non-interacting single species density of states. 

We now decouple the four fermion interac- 
tion term via a Hubbard- Stratonovich trans- 
formation in the exchange channel, to obtain 
Z = J DtpiDil>Dhexp[- J dxC{^ ,i>,h)], where 

C = V f [&t + H + h] V> + ^Tr[M+]. (11) 

Here, h is an 8 x 8 hermitian matrix, which we write as 
h = M ® Q, where M is a 2 x 2 hermitian matrix in 
sublattice space and Q is a 4 x 4 hermitian matrix in 
flavor space. The gapped states (Ref . [U-lfij] ) correspond 
to taking M — At 3 . Integrating out the fermions yields 
Z = J D(Q)D(A) exp(- / dxC[A(x),Q(x)], where 

C(A,Q) = -Trln [d t + H + Ar 3 Q] + ^Tr[Q 2 ]. (12) 

The SU(4) flavor invariance manifests itself in an exact 
SU(4) flavor degeneracy of the many body states. Upon 
minimizing the action (|1 2|) in a saddle point approxima- 
tion, we find Q 2 = 1, and A = A exp(— 2/Xi/q), where 
A w 0.4 eV is the bandwidth for the two band Hamilto- 
nian. This gives the mean field Hamiltonian, Eq.©. 

We note that instead of decoupling the interaction in 
the excitonic channel h = T3 Cg> Q, we could have chosen 
the channel h = n,2 <8> Q- This choice would lead us to 
the nematic state of Ref0, which is gapless, but breaks 
lattice rotation symmetry. However, the nematic state 
is higher in energy than the gapped states at the saddle 
point level, so we will concentrate on the gapped states, 
and specifically on the lifting of the accidental degenera- 
cies by thermal and zero-point fluctuations. 

Our symmetry analysis, involving multiplcts 
(M>,M<) for different matrices Q, could also be 
applied to the nematic stated. However, the fluctuation 
analysis cannot be perfomed because the T38Q mode has 
negative rigidity, i.e. the nematic mean field is unstable. 



VI. LIFTING THE DEGENERACY: 
ZERO-POINT FLUCTUATIONS 

We first analyze the case of zero temperature, when 
the degeneracy is lifted by zero point fluctuations. The 
most important fluctuation modes are those that are soft- 
ened under RG. In BLG, this means the 'L' modes t^SQ, 
which describe fluctuations longitudinal with respect to 
the order parameter in sublattice space, and also the 'T' 
modes T1.2SQ, which describe fluctuations transverse to 
the order parameter in sublattice space. In that, 5Q is 
an arbitrary 4x4 hermitian matrix. 

We therefore expand the action in Eq. (TT2")) to quadratic 
order in the fluctuation modes r Q ® SQ a , a = 1, 2, 3, to 
obtain 



(13) 



Here, Latin indices k = 1...N refer to fermion flavor, 
whereas Greek indices a, (3 = 1, 2, 3 refer to the Pauli ma- 
trices r a that parameterize the fluctuations in sublattice 
space. The matrix K is defined by 
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where we have introduced the polarization operator 
n?/( W ,q) = J ^^Tr^Gi^p+yrpGiip-)). (15) 

It is convenient to choose a diagonal background state 
Q = Ci^r/i where Q = il) so that the Greens function 
takes a form diagonal in the flavor space, 



Gi(p±) 
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The trace in Eq. ([T5"j) goes over sublattice indices, but not 
over flavors. 

The matrix K is positive definite, so we may integrate 
out fluctuations to obtain an expression for the fluctua- 
tion contribution to the free energy, 



F fluct = i^^lnK^( W ,k), 



(17) 



aij cjk 



where we took into account that the only contribution 
comes from the diagonal terms, a = fi, i = I, j = k. We 
now subtract the fluctuation energy of the (4, 0) QAH 
state from that of the (2, 2) state, to obtain 
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SF - i*fl uct ,(4,0) - -Pfluct,(2,2) = 4 E ln 1 T T~la a ' ( 18 ) 
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where 11^% and n?^ are defined by Eq. (fT5"|) . with 
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(Ci)O) = (1, 1) and (1,-1) respectively: 

nf>,q)=||^iTr(r Q G > ( P+ )r Q G > (p_))(19) 
D5^(a;,q)=y^i'&(r Q G>(p + )r a G<(p_))(20) 

G>(<)(e ' PH ,e-H (p)TAr3' 

where we used a shorthand notation p± = (e ± ^w, p ± 
iq). To analyze the effect of competition of different 
modes in full detail, below we compare the fluctuation 
energy for the states of different type (M>, M<). 

To evaluate the difference of fluctuation energies, given 
by Eq. (fT5)) , it is convenient to rewrite it as 



SF 



a = l 



1 | , , < ! 

11 >> 



TTOQ 
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(21) 



Below we evaluate the differences of polarization func- 
tions n"^ — n>< , and find that different modes, L and 
T, yield contributions of opposite sign. 

In particular, we find that Tl^ > — is positive, i.e. 
the L-modes favor the (2,2) state. This effect of longi- 
tudinal modes is well known in the topological insula- 
tor literature 2 -. In contrast, the differences II"" — II" < 
with a = 1, 2 are negative. Thus, the T- modes, which 
are unique to BLG, favor the (4,0) state. We evaluate 
Eq. (fT5)l . and find that the T- modes dominate the free 
energy, favoring the QAH state. 

To proceed with the analysis of the quantities ]!"> — 
II"" , it is convenient to define e± = e ± cj/2, and z± = 
|p ± iq| 2 /2m. In this compact notation, we have 



where / ... = f 



and D aa (u),q) 



II"<(u;,q). The integral over e may be performed ex- 
actly by the method of residues, to give 

diJcPqcPp ( 1 
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Du(w,q) D 22 {uJ,q) 
J- + 



^33(w,q)y u, 2 + (£++£_) 2 ' 



where £± = W z± + A 2 . The integral over p may now 

be performed with logarithmic accuracy. The dominant 
contributions come from £± » 0, and may be evaluated 
as 
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where we have used the pseudo-polar coordinates r 2 



(q /2m) and have assumed that r> A. 



We now have to calculate the various functions D aa . 
We will calculate these quantities analytically with log- 
arithmic accuracy. We begin with the definition D aa — 
j + II" <, where the polarization functions are defined 
in Eqs. ([T9|) .([20 |) . We note that the polarization func- 
tions n^ Q are logarithmically divergent at small u, small 
|q| 2 /2m and A = 0. The coefficient of the logarithm 
can be extracted by setting ui, q, A = in the integral 
in Eqs. (fTi?l) . (j2T)|) . and introducing an IR cutoff r, where 

r 2 = to 2 + (|q| 2 /2m) 2 , and we assume r > A. In this 
manner, we obtain 
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F aa (e,p) = A 2 Tr(r a T~3T a T 3 ) + ATr( T3 )(«£+) 
+ATr(T a H {p + ^q)T a T 3 ), 

where / ... = j* 7§fyr — Terms in Eq.([22]l linear in 
A must vanish, since the fluctuation energy should 
be invariant under sign changing A — > —A. Techni- 
cally, the vanishing of terms linear in A follows because 
Tr(r a T 3 T a ) — 0, and Tr(r a i?o T a''3) = 0. As a result, the 
first term in F aa (e, p) (at order A 2 ) is the only term that 
survives. We can substitute the expression in Eq.(| 



into Eq. (fl"8|) and expand the logarithm in small A 2 , to 
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(23) 
(24) 
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where we have introduced the notation z 



(26) 
(27) 

and 

p| e »«p _ p x _|_ ^he integrals may be straightfor- 
wardly performed by changing to the pseudopolar coor- 
dinates (p,ip,9 p ), where e = pcosip, z = psinip, and 8 P 
was defined above. The integral goes over < 9 p < 2tt, 
< ip < 7r and r < p < A. Integrating in turn over 9 P , ip 
and p, we find 



JPT 

2m 



Dn = D 2 



A 4 n r ' 



1 A 
D 33 = ---ln-. (28) 



We now recall the relation A 1 = \ Vq In A/ A (the gap 
equation), and substitute it into Eq. ([23|) . to obtain 
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This integral can be evaluated using the substitution x 
m 3, giving 



into account L-modes only, is given by a sum over Mat- 
subara frequencies, 



SF = 8 A 2 up r^'r'Yf'- 3 ^. PO) 
u io ln(A/A) + a; v 1 

Evaluating the integral, we obtain a negative value 
SF = 8(-3 + 41n2)A 5 V ln-^ w -1.82A 2 j/ In (31) 

which favors the QAH state. 

It should be noted that the difference in energies be- 
tween the (4,0) and (2,2) manifolds is of the same order 
as the mean field energy, so the mean field plus fluctu- 
ations analysis is ill controlled. However, it provides us 
with an intuition about the splitting between manifolds 
of different signatures, and we believe the qualitative de- 
tails of the fluctuation splitting are reproduced correctly 
by this analysis. 

We note that our fluctuation analysis included only 
those modes that correspond to weak coupling instability 
in BLG. We could also have included Stoner modes in our 
fluctuation analysis. These would produce an additional 
contribution 



W^Stoncr — 8A 2 i/ 
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(32) 
(33) 



where H^ < is defined by Eq. (|2"0f with a 

TOO 



P = 0, i.e. 

with r a = Tp = 1. Now, since IKJL is not log diver- 
gent, we can take D 00 = 1/A with logarithmic accuracy. 
We then obtain a contribution dFs% oaeT = 4A 2 j/olnA/A, 
which is sufficiently large to change the sign of the result 
Eq . (f3Tj) . However, this calculation, which neglects corre- 
lation effects, is likely to strongly overestimate the effect 
of Stoner modes, and therefore we believe that Stoner 
modes should be left out of the fluctuation analysis. 



F &uct = \t £ ln (j + U i! k> k ) 



(34) 



where uj n — 2imT. 

We will perform a long wavelength expansion of 
Ilf?(a;,k). At zeroth order, we note that at uj, k = 

>< " iivi ix >> 



the values of H^ < and are given by 



H^k^O) = \ J |^tr (r 3G> r 3G> ) 



(2tt) 3 

e 2 + z 2 - A 2 d 2 pde 

( £ 2 + Z 2 + A2)2(2^)3' 



nf < (^ = 0)= 1 -J |^tr (T3G>r 3 G< 



d 2 pde 



(£ 2 + 2 2 + A 2)2 ( 2?r )3 ' 



where G>«) = l/(ie — H (p) =p At 3 ). To distinguish 
Goldstone modes from gapped modes, it is convenient to 
recall the gap equation 



1 



d 2 pde 



A 2 (2tt; 



3 ' 



(35) 



Hence, we have j- + n? > 3 < (0) = 0, which corresponds to a 
Goldstone mode, whereas in the case of we have 



i + n 3 > 3 > (0) = / 



2A 2 



d 2 pde 



{e 2 + z 2 + A 2 ) 2 (2tt) 3 ' 



(36) 



which is manifestly positive. Thus, Goldstone modes ex- 
ist only in states (Af>, M<), where M> ^ and M< ^ 0. 

The free energy, Eq. (|34|) . evaluated at leading order in 
a long wavelength expansion around w, k = 0, is given by 
a sum 



VII. LIFTING THE DEGENERACY: THERMAL 
FLUCTUATIONS 

Thermal fluctuations are dominated by gapless Gold- 
stone modes, which are present only in the states that 
break SU(4) symmetry. In a state (M>, M<), there are 
M > M < Goldstone modes. Thermal fluctuations due to 
Goldstone modes allow a state to gain entropy, and since 
the (2,2) states have the most Goldstone modes, they 
have the highest entropy. It may thus be expected that 
the (2,2) states dominate at sufficiently high tempera- 
ture. 

Below we present an analysis showing that this expec- 
tation is correct. Since gapless fluctuation modes appear 
only in the L-mode channel Sh oc t 3 SQ, it is sufficient 
to restrict our attention to the L-modes. The general ex- 
pression for the fluctuation part of the free energy, taking 



Ffluct = TJ2 M > M < Haul + bk 2 ) (37) 

+ i(Af 2 + Ml) Ha'ul + b'k 2 + c), (38) 

where the first term is the contribution of the gapless 
modes (originating from n?^), while the second term is 
the contribution of the gapped modes (originating from 
n?^). The coefficients a, a' , b, b' are obtained by Taylor 
expanding n 33 (w,k) in small u and k, while c is given 
by Eq.©. 

To simplify the sum over Matsubara frequencies, it is 
convenient to define the quantity f(u) = T ^ w k ln(cj 2 + 
u 2 ). We can evaluate f(u) by first taking the derivative 



df x 1 

-r = T Y - — 



oj n ,k 



+ c.c. = coth — , 

ILOn + U 2T 



7 



and then integrating it over u to obtain 

f(u) = 2Tlnsinh^ = 2Tln (l - e - u/T ) +u-(21n2)T. 

Plugging this identity into the sum (I57|) . we see that the 
contribution of the gapped m odes is exponentially small 
at low temperatures, T <C \J c/a' ~ A, while the sum 
over gapless modes gives a negative contribution of a 
power law form, 

Ffl uc t = 2M>M<Tln (l - e -«W/ r ) , (39) 

k 

where u = 6/a ~ y A/to. Evaluating the integral, we 
obtain an estimate 

F fluct ~ -M > M<( i / /A)T 3 , (40) 

which describes the free energy gain due to thermal fluc- 
tuations of Goldstone modes. 

We see that the gapless Goldstone modes dominate 
the finite-temperature fluctuation contribution to the free 
energy. These modes lower the free energy (by increasing 
entropy). Since the number of gapless modes M > M < is 
maximal for the (2,2) states, these states are entropically 
favored by thermal fluctuations. 

What is the outcome of competition between the zero- 
point fluctuations and thermal fluctuations? In Sec I VII 
we found that at zero temperature the (0,4) QAH state 
is energetically favored by zero point fluctuations of the 
modes "softened" under RG. At the same time, the zero- 
point fluctuations of other modes, such as the Stoner 
modes, may have an opposite effect, favoring the (2,2) 
state. In the event the zero-point fluctuation energy is 
dominated by such non-soft modes, the (2,2) state will 
be realized in the entire temperature interval where the 
system is unstable to gap formation. 

A more interesting situation may arise if the zero- 
point fluctuation energy is dominated by the RG-softened 
modes, favoring the QAH state at zero temperature. In 
this case, given the opposite effect of zero-point and ther- 
mal fluctuations, we have to consider the competition 
between the QAH and (2,2) states. Since the thermal 
fluctuation energy (jlO]) vanishes at T = 0, we expect 
that zero point fluctuations will dominate below a cer- 
tain temperature T», above which thermal fluctuations 
will dominate. If T* < T c , where T c » A(T = 0) is 
the critical temperature for gap opening, then a QAH 
state will be realized at low temperatures < T < T», 
whereas a (2,2) gapped state will be realized in the in- 
terval T* < T < T c . In contrast, if T* > T c , then the 
QAH state will transition directly to an ungapped state 
at T = T c via a second order phase transition, and the 
(2,2) state will not be realized. 

A rough estimate of the temperature T* can be ob- 
tained by comparing the free energies (j4*U|) and ([3Tj) . 

T 3 A 
£Ffluct,(2,2) 4 ^0"^"' ^auct,(o,4) 1.82^ A 2 ln— , 

(41) 



indicating that the scale for T» is comparable to the tem- 
perature T c at which the gapped state forms. A more de- 
tailed analysis of temperature-driven transition between 
the QAH state and (2,2) state is beyond the scope of this 
work. 



VIII. EXPERIMENTAL SIGNATURES OF THE 
QAH STATE 

We now discuss experimental tests of the QAH state. 
The clearest experimental signature would be detection 
of the quantum Hall effect at zero external magnetic field. 
However, detection of this effect requires four probe mea- 
surements performed on a sample of BLG that is suffi- 
ciently clean and at sufficiently low temperatures as to 
exhibit spontaneous gap opening 5 . Such measurements 
have not yet been performed. Moreover, detection of this 
effect could be complicated by the formation of domains 
with opposite signs of A. Different domains will have 
opposite a xy , so the Hall conductance of a macroscopic 
sample will average to a value near zero. However, if 
there is percolation of edges, there will be a non- vanishing 
two-terminal conductance of order e 2 /h. 

Alternative experimental tests of the QAH state may 
be performed by examining the electronic compressibility 
in weak magnetic fields. When the chemical potential sits 
near the missing Landau level in Fig.([T]), there should be 
a gap that extrapolates to a non-zero value as B — >• 0. 
This effect will be seen at either y = 4 or v = - 4 if there 
is only one domain, and at v = ±4 if there are multiple 
domains. 

The gap at v — 4 will be strengthened by the mecha- 
nism outlined around Eq.(|5]), however, a signal at v = —4 
will be seen only if the QAH state is intrinsic, rather than 
field induced. An incompressible region at v = —4 com- 
bined with a gapped state at B = can thus be taken as 
a diagnostic for a QAH state at B — 0. The filling factors 
v = ±4 are not equivalent because the QAH state breaks 
particle-hole symmetry in magnetic field. 

Another experimental signature is a phase transition 
at filling factor v = and finite B from a QAH state 
to the Quantum Hall Ferromagnet (QHF) states that 
are expected to form at large magnetic fields^ 3 .. Such a 
phase transition would not be seen if the dominant state 
at small B was of (2,2) type, since the (2,2) states are 
smoothly connected to the QHF state. 

An incompressible region at v = ±4 that occurs at 
anomalously low magnetic fields, such that the features in 
compressibility at other integer v values are washed out, 
was found in recent experiments that employed a capaci- 
tance scanning probe to study suspended BLG samples 9 . 
In transport measurements^ performed on the same sys- 
tem, a state with finite two-terminal conductance of or- 
der e 2 /h was found at zero field, which at a finite B 
field undergoes a transition to an insulating state. These 
measurements are all compatible with the QAH state, 
however, since there is as yet no four-terminal measure- 
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ment, it is not possible to say for certain whether a QAH 
state has been observed. 

In summary, our symmetry classification of the vari- 
ous gapped states proposed for BLG singles out the QAH 
state as the only gapped state not breaking any contin- 
uous symmetry. We have investigated the fluctuation- 
induced splitting of the gapped states, and concluded 
that at zero temperature and zero field, the leading in- 
stability is to the QAH state. We have discussed the 



phenomenology and experimental signatures of this state, 
and have shown that it can be stabilized by weak external 
magnetic field. 
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